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Considering the long range Coulomb interactions between large polarons in dielectrics, we propose 
a model for their crystallization when no bipolarons are formed. As the density increases, the 
melting is examined at T = OK. One possibility is the delocalization towards a liquid state of 
polarons. However, we show that this cannot happen if the electron-phonon coupling is larger than 
some critical value. The other competing mechanism is the dissociation of the polarons themselves, 
favored owing to their large mass at strong coupling. Finally, we propose a phase diagram for the 
insulator-to metal transition as a function of the density and electron-phonon coupling. 



Any theory of band structure for interacting fermions is 
based on the fact that the long range part of the Coulomb 
repulsion is screened and contained in collective modes, 
called plasmons [Q. However, for an idealized low density 
system the long range Coulomb repulsion dominates the 
kinetic energy, and Wigner Q has shown that the carriers 
crystallize in that case. This does not occur for electrons 
in usual metals, owing to the high density of carriers. 

For doped insulators, the carriers' density induced by 
doping can be extremely low. These materials can be 
classified into two types: polar and non-polar. Non-polar 
compounds are the usual covalent semiconductors which 
are easily doped up to the insulator-to-metal transition. 
Such materials are characterized by their dielectric con- 
stant K, which does not depend on the frequency (the 
static dielectric constant Eg and the high frequency di- 
electric constant £oo are close together k « Eg ~ £oo)- 
The insulator-to-metal transition is due to the screen- 
ing and ionization of the localized impurity states (the 
effects of disorder are neglected here). Mott |^ has 
given a criterium to evaluate the critical density above 
which the ground-state is metallic. Given the radius 
Rimp. = h^nje^m* of the impurity bound state (m* is 
the effective mass of the carrier), the system is metal- 
lic for doping densities n such that n > , where 
{n^^y^^Rimp. ~ 0.25. Mott assumes that the long range 
Coulomb interactions between liberated carriers are also 
well screened as soon as n > ric. This is correct because 
at these densities the kinetic energy of the carriers is 
much higher than the Coulomb energy 1/kRs {Rs is re- 
lated to the density by 1/n = AttRI/S), which is reduced 
by K (which can be large), so that we are sure that the 
Wigner localization would occur only at densities much 
lower than the critical Mott density. 

For polar insulators, such as some transition metal 
oxides, the situation is more complicated owing to the 
formation of polarons. When the static (Es) and high 
frequency (Eoo) dielectric constants have an appreciable 
difference of magnitude, a free carrier is bound to the po- 



larization that it creates around itself. The couple made 
of the carrier with its surrounding lattice polarization 
is called a polaron For the single polaron prob- 

lem, one defines the dimensionless electron-phonon cou- 

1/2 

pling a — (^m* /2h^LJir) e^/i. The constant e, which 
is responsible for the polaron formation, is defined as 
l/i = I/Eoo — l/ssj and ujir is the phonon frequency 
of the optical longitudinal mode. The polaron problem 
has been extensively addressed in many theroretical stud- 
ies However, Feynman Q has proposed a variational 
method using path integral calculations giving with a 
very good accuracy the ground state properties for any 
value of a. The idea is to use a simplified variational 
Hamiltonian to calculate quantities of interest: 
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Here the electron with coordinate x and momentum p, 
is coupled to a fictitious particle with coordinate X and 
momentum P which simulates the polarization field. K 
and M are the variational parameters. As pointed out by 
Feynman, eq.(^ well imitates the polaron. In particular, 
M gives a good estimate to the true effective mass of the 
polaron : Mp « m* -\-M ^ and the polaron radius can also 

be defined as Rp « \J n/2Mp{K/Mpf''^ . Table 1 gives 
the rest mass Mp obtained by Feynman and Schultz 0| 
for different values of a. We notice that for large a, the 
polaron mass is much larger than the mass of the free 
carrier m*. 

If one considers now two polarons in the system, each 
electron (or hole) is attracted by the polarization in- 
duced by the other. On the other hand, there is also 
a direct repulsion between them. The resulting competi- 
tion yields two different situations depending on the ratio 
"H — Soo/^s- Many theoretical calculations |^ have shown 
that if 77 < « 0.1, a bound state between two polarons 
is formed, which is called a dielectric bipolaron. If this 
condition is not satisfied, the Coulomb repulsion domi- 
nates and both polarons repel as l/e^c?, where d is their 
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relative distance. We will assume from now on that we 
are in this last situation, which is often verified in real 
materials. All along this letter, we will consider materials 
which verify 0.1 < rj < 0.5. 

On chemically doping at finite density polar dielectrics 
with 77 > 0.1, we are faced with three intricated prob- 
lems: (i) the polaron formation, (ii) their binding to the 
impurity atoms (known to be the F-Center problem [^), 
and (iii) the Coulomb repulsion between polarons if they 
are mobile. As already pointed out by one of us , the 
last point (iii) is essential to understand the physics of 
such systems. Since the polaron mass can be large (Table 
1), the kinetic energy of a set of polarons which behaves 
as 1/MpR^ can remain small compared to the Coulomb 
energy l/ssRs even at high density and high dielectric 
constant. To simplify the problem and focus mainly on 
the effects of the long range Coulomb interactions, we 
shall concentrate in this letter on an idealized system 
of polarons at finite density and interacting in a jellium 
of equal density (and opposite charge). At low density, 
they crystallize jl^l as the electrons in the Wigner crys- 
tal. We will now address the problem of the melting of 
such a crystal as the density increases. 

The strong electron-phonon coupling limit a — > cx), 
or equivalently the low phonon-frequency limit ujir — > 
can be first examined. It is well-known that in this 
limit, the polaron is formed by an electron localized in 
a quasi-static polarization potential-well, which behaves 
as —e^/ir at large distance, and is almost constant at 
short distance. The mass Mp of the polaron becomes in- 
finite as ujir vanishes, while its energy Ep and radius Rp, 
which are finite, only depend on the dielectric constants 
and effective mass m* of the electron. When the density 
increases, the polarons are unable to delocalize owing to 
their infinite mass. Thus, at some critical density, the 
electronic states localized within the polarization-wells 
must disappear for the same reason as the localized impu- 
rity states at the usual Mott transition: they are screened 
by liberated electrons. Making the same assumptions as 
Mott in the case of impurity states, we propose the fol- 
lowing qualitative criterium for the dissociation of the 
polarons, valid for LUir ~ 0, i.e a —^ 00: 
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Rp ~ Const. 



(2) 



The ratio e/Soo takes into account the different dielectric 
constants appearing respectively in the calculation of the 
Thomas- Fermi length (Eoo) and in the polaron radius (e). 
The resulting state occuring for n > will be discussed 
later on. Numerically, on taking the values Eg = 30, 
£00 — 5, m* = nie and Const = 0.25 [pl| , one obtains 

TT-c ~ lO^^CTO"^. 

The above consideration of the insulator-to-metal tran- 
sition based on the polaron dissociation is a priori no 
more valid at finite electron-phonon coupling, since in 



that case the polarons could delocalize to form a liquid 
state instead of being dissociated. These two compet- 
ing mechanisms for the insulator-to-metal transition are 
now more precisely examined at finite electron-phonon 
coupling and zero temperature. To modelize the situa- 
tion at finite polaron density, we add to the host material 
a density n of electrons plus a rigid compensating jellium 
of same density. To account for the strong Coulomb in- 
teractions in a kind of mean-field approximation, start- 
ing from the low density limit, we take as a simplified 
model, a single electron localized within a uniformly pos- 
itively charged and polarizable sphere of radius Rs, which 
is nothing but the Wigner crystal of polarons (see ref. jl^] 
for details). As done by Wigner for the electron crystal, 
we neglect the dipole-dipole interactions between adja- 
cent spheres. The hamiltonian describing such a system 
is then given by: 
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where D+ and D~ are the electric displacement due 
to the jellium and to electron respectively. Pir is the 
polarization. We now introduce the polarization part 
which only responds to the electron motion : P,:r = 



(l/47re)D+ and P, 



and easily obtain: 
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Co = —9e'^ /IOEsRs is a constant part of the energy at 
fixed density, and UJ^/ = /m*£sRl is related to the 
usual plasmon frequency = Anne'^ /m* by the relation 
= Wp/SSs. Equation (||) is nothing but the hamilto- 
nian of a polaron, the electron of which is bound in a har- 
monic potential. Taking the limit i?s — s- cx), one recovers 
the usual single polaron hamiltonian, with Cq ~ luw — 0; 
D+ = 0, and Pir = P.,r- 

As previously proposed by Nozieres et al. ||l^,|l^, the 
Lindemann criterium allows a good qualitative discussion 
for the melting of the electron Wigner crystal. The essen- 
tial idea is that the crystallized state is unstable when the 
quantum fluctuations (zero point motion) of the electron 

positions exceed some critical value. More precisely, the 

1/2 

crystal melts when <^(5a;^) /Rg > 5, where 5 is a phe- 
nomenological constant usually taken as (5 « 1/4. For 

1 /2 

the Wiffner crystal, {5x^) = n/2m* {e^/m*Rl) '\ and 
the critical value of Rg below which the crystal melts |l^] 
is [H : i?^ = 64ao (ao is the Bohr radius). 
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The apphcation of this criterium to the polaron crystal 
is more complicated owing to the composite nature of the 
polaron. To be able to delocalize the polarons towards 
a liquid state, the quantum fluctuations for each single 
polaron moving as a whole, i.e. the motion of the elec- 
tron together with the surrounding polarization, must 
become large in comparison with Rs. On the contrary, 
if the relative fluctuations of the electron with respect to 
the polarization become large, the polaron dissociation 
is favored. The Feynman model yields a natural way to 
evaluate the behavior of such quantities: since the po- 
larization field is replaced by a particle with coordinate 
X, one easily calculates both the quantum fluctuations 
of the center of mass R = {m*x + MX)/{'m* + M) and 
the relative coordinate r = x — X . Thus, we have now 
two different Lindemann criteria: 

1/2 

(i) {5R^) /Rs > 1/4: which allows the melting of the 
crystal towards a liquid state of polarons ; 

1/2 

(ii) ((5r^) /Rs > 1/4 for the polaron dissociation. 
Both quantities are calculated by the same variational 

treatment, just taking instead of equation (2), the new 
Hamiltonian: — Hq +7x^/2, where 7, K and M are 
the new adjustable parameters. Once the lattice degrees 
of freedom are integrated out , the action for the model 
(0) is given by: 



S^^J dtx^(t) + J dtx^t) 

~V8J J R^F^' 



(5) 



where the integrals are calculated over a time interval T, 
h = 1 and the frequency unit The trial (quadratic) 

action corresponding to is now 



dtx^(t) + - / dtx\t) 



C 



dt ds [x{t) — x{s)] e 



(6) 



where — K/M and C ~ Kw/4. 

At low density, the two degrees of freedom r and R are 
uncoupled, and their frequencies are respectively LOint and 
ujext- These are sketched in Fig. Q as functions of n. The 
two corresponding quantum fluctuations are then given 
by 5r w {h/2m*ujintf''^ and 5R « {n/2Mpujext)Y'^ ■ 
Fig.^. (a) shows the numerical values of the ratios (i) and 
(ii). One observes that for fixed a, while {Sr"^)^/^ /Rs in- 
creases with increasing density, {SR"^)^^^ / Rg is bounded 
from above. In addition, its maximum value decreases 
with increasing a, so that above a certain critical value 
a* , the value 1 /4 is never reached and the crystal melting 
is prevented. 

The situation encountered at infinite electron-phonon 
coupling (eq.(2)) is thus continued down to finite values 
of a. For a > a*, the liquid state of polarons cannot 



exist. Physically, u>ext. being limited by the frequency of 
the phonons themselves {iVext. < i^ir, see Fig. |^), our 
result means that, when the electron-phonon coupling is 
too large, the polarization cannot follow the increase in 
kinetic energy induced by the doping through tow- The 
quantum fluctuations are then transferred to the internal 
degree of freedom. For that reason, the polarons dissoci- 
ate at the transition. 

Let us focus more precisely on the shape of 
{SR'^y^^/Rs in Fig. |(a). Starting at low density 
{Rs 00, uJw ~ 0), the extra term m*Cij'^x1/2 is a 
perturbation with respect to the single-polaron energy, 
and the polaron properties are almost unchanged as long 
as its radius is much smaller than the average distance 
between carriers. To a first approximation, we can then 
neglect the composite nature of the polaron and treat it 
as a rigid particle of mass Mp (the single polaron mass 
from Table 1) moving in a weak harmonic potential with 
a frequency Lu^xt — 

{m*/Mpf/^. The average quan- 
tum fluctuations in the jellium potential are given by the 
following expression 
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which fits very well the low density part of {SR"^)^^^ /Rs 
(right side of Fig. ||.(a)). This gives the critical cor- 
responding to the first Lindemann criterium (i) in terms 
of R'^ = 64ao 



Ri'^ = Rl 



{Mp/m* 



(8) 



If the mass of the carriers is strong enough, despite the 
static screening of the charges, the crystallisation is fa- 
vored in comparison with the ordinary Wigner lattice. 
On the other hand, (Jr^)^/^ in (ii) is nothing but the 
polaron radius, so that (ii) gives : i?c"^ = ARp. 

At higher densities, the frequency of vibration of the 
polaron uJext = towiTn* /Mp)^/"^ increases, due to the 
tightening of the jellium potential. Thus, if the den- 
sity is such that uJext ~ i^ir, i-e. cjw ~ i^iriMp/m*Y^'^, 
the quenching of the polaron frequency is translated into 
an enhancement of the polaron mass, which can be de- 
scribed through a function Mp{Rs). This is observed 
numerically, as shown in Fig. |[(b) (the rest mass is 
recovered for Rs 00). As a consequence, at high 
density, the quantum fluctuations of the polaron will 
be {n/2Mp{Rs)ujirY''^ for the left part in Fig. |.(a), 
while eq.(|^) is valid for the right part. This gives rise 
to the bell shape of {SR"^)^^^ / Rs- This allows a simple 
evaluation of a*: the maximum value of {SRY^^/Rs 
will be obtained approximately by evaluating the low 
density formula eq.(^ at ujw ~ LUir{Mp/m*y^^ . If 
we equate the result to the melting value 1/4 we get 
{Mp/m*y/^a « 16(es/e). As an example, taking eg = 30 
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and = 5, and using the values for Mp of Table 1, one 
obtains a* ~ 9. The variational result, without such 
approximations, gives a* ~ 7.5. 

The phase diagram Fig. || gives the calculated critical 
radii for crystal melting (i) and polaron dissociation (ii) 
as functions of the electron-phonon coupling a. Though 
the transition radii decrease linearly with Es and e, the 
general shape of the diagram does not change in the stud- 
ied range 0.1 < Sao/ss < 0.5, and so does the critical 
value a* {a* « 6 — 9). 

For small a the transition corresponds to the melting 
of the usual Wigner crystal of electrons weakly renormal- 
ized by the electron-phonon coupling. The corresponding 
metallic region ((I) in Fig. 2) has been previously exam- 
ined by Mahan [|l5| . Extension of his calculations in the 
metallic state have been performed by Devreese et al. at 
intermediate coupling, and more recently by ladonisi et 
al. at strong coupling |l5| . 

For intermediate coupling (say 5 < a < a*), it can 
be expected that the polarons remain well defined, yield- 
ing a liquid state of polarons below the transition (re- 
gion (II) in Fig. 2). The physics in that case may be well 
described by the Holstein model for non-interacting po- 
larons 

For strong coupling (a > a*), the crystal cannot melt 
towards a polaron liquid. At the transition, the po- 
larons dissociate and the resulting metallic state (region 
(III) in Fig. 2)) will present both strong electron-phonon 
and electron-electron interactions. A theory for such a 
system has been recently proposed by Castellani et al. 
p8[ : it presents phase separation due to the electron- 
phonon coupling, which is prevented by the Coulomb 
repulsion. The resulting effects are strong charge den- 
sity wave fluctuations, and a superconducting instabil- 
i ty. Another possibility, already proposed by one of us 
|lO| is that only part of the polarons dissociate to quasi- 
free electrons at the critical density. In that case, just 
above the insulator-to-metal transition the system would 
present a kind of electronic pseudo-separation between 
localized polarons and quasi-free electrons. However, at 
higher density, this pseudo-separation would disappear 
owing to screening effects, and in that case the system 
would be well represented by the theory of Castellani 
al. [p^ . This possibility will be examined in a further 
publication. 

According to our theory, the possibility of polaron dis- 
sociation at the insulator-to-metal transition in real ma- 
terials depends on the value of a, which is generally lower 
than a* « 6 — 9. However, any mechanism which in- 
creases the polaron mass Mp compared to the present 
3D case, will shift a* to lower values and make the dis- 
sociation of polarons physically reasonable. Moreover, 
the critical density for this dissociation will be increased 
compared to the present 3D case (n^^ ~ IO^'^^cto"^, see 
eq. (2)). Two examples of such mechanisms can be men- 
tioned: (a) a strong anisotropy of the effective mass m* , 



which corresponds to considering two dimensional po- 
larons; and (b) the formation of a polaron in a Mott- 
Hubbard insulator, which is associated to a spin-polaron 
]l9| , which mass increases exponentially with its size. 

In conclusion, we have proposed a model which takes 
into account in a self-consistent way the effects of the 
Coulomb interaction for a system of polarons at finite 
densities. A new mechanism for the insulator-to-metal 
transition has been put in light, which corresponds to the 
polaron dissociation at strong electron-phonon coupling. 
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Mp/m* 


1.8 


3.6 


13.2 


59.2 
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TABLE I. Polaron effective mass as given by Feynman. 



FIG. 1. The two frequencies LOint. and uiext. as functions of 
the density are represented in arbitrary units for a = 10. The 
external frequency corresponding to the quantum fluctuations 
of the polaron as a whole saturates and is necessarily lower 
than ujir- The critical density ric for the polaron dissociation 
(see text) is also quoted. 

FIG. 2. (a) The ratios {SR^y^'^/Rs (bold line) and 
{Sr^y^^/Rs (dotted line) versus Rs/cio at a = a* (see text), 
for too = 5, = 30 and rn* = nie. (b) The corresponding 
effective polaron mass Mp in units of m* . 

FIG. 3. Phase diagram for e^o = 5, Es = 30 and m* = me- 
Crystal melting (bold line) and polaron dissociation (broken 
line). The dotted line is the critical radius given by eq. (8). 
The arrow marks the result for the modified Mott criterium 
eq. (2). For regions and III: see text. 



5 



0.6 
0.5 



0.4 



0.3 
0.2 



0.1 







1 



(a) 



10000 



10 100 1000 10^ 

s 




100 1000 10^ 



s 



104 



103 



102 



101 



1 



Wigner 






Crystal 








r oiaron 




^» 
% * 


Crystal 




% \ 
% \ 
% % 
% « 

II \'n 

■ 
1 
1 






1 \ 


III 




Metallic State 











10 
a 



15 



20 



